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Abstract. We prove regularity and stochastic homogenization results for cer- 



tain degenerate elliptic equations in nondivcrgence form. The equation is re- 



quired to be strictly elliptic, but the ellipticity may oscillate on the microscopic 
scale and is only assumed to have a finite dth moment, where d is the dimen- 
sion. In the general stationary-ergodic framework, we show that the equation 
homogenizes to a deterministic, uniformly elliptic equation, and we obtain an 
explicit estimate of the effective ellipticity which is new even in the uniformly 
elliptic context. Showing that such an equation behaves like a uniformly el- 
liptic equation requires a novel reworking of the regularity theory. We prove 
deterministic estimates depending on averaged quantities involving the distri- 
bution of the ellipticity, which are controlled in the macroscopic limit by the 
ergodic theorem. We show that the moment condition is sharp by giving an 
explicit example of an equation whose ellipticity has a finite pth moment, for 
every p < d, but for which regularity and homogenization break down. While 
these results are new even for linear equations, we prove them in the fully 
nonlinear context. In probabilistic terms, our homogenization results corre- 
spond to quenched invariance principles for controlled diffusion processes in 
random media, including linear diffusions as well as diffusions controlled by 
one controller or two competing players. 



1. Introduction 

We prove stochastic homogenization and regularity estimates for fully nonlin- 
ear elliptic equations in nondivergence form without the assumption of uniform 
ellipticity. The equations we consider are strictly elliptic but may have elliptici- 
ties which are arbitrarily large and oscillating on the microscopic scale. We derive 
new (deterministic) regularity estimates and show that, under the assumption that 
the dth moment of the ellipticity is finite, such a degenerate equation homogenizes 
in the macroscopic limit to an effective equation which is uniformly elliptic. Our 
analysis yields an explicit estimate for the effective ellipticity which is new, to our 
knowledge, even in the linear, uniformly elliptic setting. In terms of probability, the 
main homogenization result, in the linear case, is equivalent to a quenched invari- 
ance principle for a diffusion in a random environment. For a nonlinear, positively 
homogeneous equation, the homogenization result gives similar information about 
the quenched behavior of controlled diffusions in random environments. 
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The simplest example of interest is the linear equation 

d 

(1.1) X "<•(-•-•)"< ' inf/CR^, d>l. 

The coefficient matrix (a^), which depends on the random parameter w, called 
the environment, is assumed to be stationary-ergodic and to satisfy the ellipticity 
condition 

d 

(1.2) A(w) |£| 2 < a ij(°^)^ < A \Z\ 2 for all £ £ K d , 

where A > is a given constant and A = X(lu) is a nonnegative random variable. 



Papanicolaou and Varadhan 25, 26] proved that, if the equation is uniformly el- 
liptic, i.e., A(w) > Ao > 0, then in the almost sure asymptotic limit e — > the 
equation (jl.ip is governed by an effective equation of the form 

where the coefficient matrix (ay) is uniformly elliptic with ellipticity A/Ao, that is, 
for every £ G M. d , 

d 

(1-3) AolCI 2 < £ < A|£| 2 - 

In this paper we extend this result to the case that A > and inf A = under the 
assumption that E[A~ d ] < oo. Due to the presence of small pockets of arbitrarily 
large ellipticity which become dense for small e, it is by no means clear at first 
glance that (jl.ll) should behave, in the macroscopic limit, like a uniformly elliptic 
equation. Our result demonstrates that this finite moment condition on A -1 ensures 
that these tiny regions of very large ellipticity may be effectively controlled and 
that (jl.ljl becomes a uniformly elliptic equation in the limit e — > 0. Furthermore, 
we show that the moment condition is sharp by exhibiting an explicit example in 
arbitrary dimension d > 1 of a linear equation, in a random environment with a 
finite range of dependence, for which homogenization breaks down but the ellipticity 
has finite pth moment for every p < d. 

Previously, in this linear setting, a result similar to our homogenization result 
has been proved in a recent paper of Guo and Zeitouni 19] using probabilistic 
methods. They give a quenched invariance principle for random walks in a random 
environment, which has an equivalent formulation in terms of stochastic homoge- 
nization of a discrete equation (on the lattice Z rf rather than R d ). They require a 
slightly stronger moment condition, namely that A -1 have a finite pth moment for 
some p > d. That homogenization occurs in the critical case p = d is new here, as 
well as the fact that it fails for every p < d. 

While our results are therefore of interest in the linear setting, we analyze much 
more general fully nonlinear equations of the form 

f(d 2 u^-,uj\ =0. 



(1.4) F(W,J,c) 
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Such nonlinear equations are more difficult to analyze than (ll.ip . due to the ab- 
sence of invariant measures- the key tool used in [2|| [2(| (and [l9l ]) to overcome the 
problem of the "lack of compactness." The results here are the first for such equa- 
tions without a uniform ellipticity assumption. Previously, the homogenization of 
fully nonlinear equations in stationary-ergodic media was proved in the uniformly 
elliptic case by Caffarelli, Souganidis and Wang Q . They introduced a new method 
based on the obstacle problem, a strategy which we also use in this paper. 

The problem one encounters when trying to homogenize (11.41) outside of the 
uniformly elliptic regime is that most of the regularity theory needed to implement 
the method of [§] is destroyed by (even tiny) regions of high ellipticity. It therefore 
seems hopeless, at first glance, to implement the techniques of [8j, since they make 
heavy use of the regularity tools. 

To overcome this difficulty prove new (deterministic) regularity estimates in 
which the dependence on a uniform upper bound for the ellipticity is replaced 
by that of its L d -norm. In particular, we prove a decay of oscillation lemma at 
unit scale. This result, and the new arguments we introduce to obtain it, are of 
independent interest. Indeed, in sharp contrast to the situation for divergence form 
equations, there are very few results in the literature for equations in nondivergence 
form which provide estimates of solutions in terms of averaged quantities. 

The estimates refine the classical regularity theory [6] and require several new 
ideas. One of the basic techniques involves using the area formula to estimate the 
size of certain "contact sets" between supersolutions and certain families of smooth 
test functions. This method is a generalization of the the classical ABP inequality 
and was previously used by Savin in his proof of De Giorgi's conjecture [28| and his 
beautiful proof of the Harnack inequality in [53] ■ In those works, supersolutions are 
touched from below not only by planes, but by translations of balls and paraboloids. 
One of the key arguments in the present paper involves touching from below with 
translations of the singular function |:e|~ q , for suitably large a > 0. 

Since the L d norm of the ellipticity is controlled on the macroscopic scale al- 
most surely by the ergodic theorem, the regularity results provide effective control 
on the solutions of (|1.4[) for small e. This allows us to homogenize the equation 
by suitably adapting the arguments of §. We expect this two-step approach to 
homogenization, in which one obtains "effective" regularity and then uses this to 
homogenize the equation, to be useful in other situations. 

As a byproduct of our analysis, we obtain an estimate for the effective ellipticity 
which is new even in the uniformly elliptic case, which states that (|1.3|) holds for 
a Ao > which, in addition to A and d, depends only on E[A _d ]. It is, to our 
knowledge, the first such bound for the homogenized coefficients of nondivergence 
form equations which is nontrivial in the sense that it is given in terms of on the 
averaged microscopic behavior of the equation rather than its uniform properties. 

As mentioned above, the homogenization result has an equivalent probabilistic 
formulation, at least in the linear quenched invariance principle for the 

corresponding diffusion in the random environment. It also provides information 
regarding the recurrence or transience of the diffusion (see [19|). If F is nonlinear 
but positively homogeneous, the fully nonlinear equation (|1.4|) is a Bellman-Isaacs 
equation which arises in the theory of stochastic optimal control and two-player 
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stochastic differential games, and the homogenization result yields similar informa- 
tion about these more general diffusion processes. Although we do not explore this 
point here, we remark that the recurrence verses transience of such controlled diffu- 
sion processes in an isotropic environment was characterized in [3|, and this result 
applied to the effective operator F, together with its proof, gives information about 
the corresponding questions for controlled diffusions in random environments. 

We now give the precise statement of our results, beginning with the modeling 
assumptions. 

The model. We work in Euclidean space R d in dimension d > 1. The random 
environment is modeled by a probability space (fi, P) endowed with an ergodic 
group r = gffid of ^-measurable, P-preserving transformations on tt. That is, 
the action r of M. d on £1 satisfies 

(1.5) F[A] = F[r y A] for every y € R d , A e 9 
and, for every A e 7, 

(1.6) r y A = A for every y e R d implies that F[A] = or F[A] = 1. 

The nonlinear elliptic operator is a map F : § x R d x D, -> R (here S d denotes 
the space of d-by-d symmetric matrices) which satisfies each of the following four 
conditions: 

(Fl) Stationarity: for every M e S d , y, z 6 R d and u £ Q, 

F(M,y,T z ij) = F(M,y + z,u). 

(F2) Local uniform ellipticity: there exists a constant A > 1 and a nonnegative 
random variable A : il — > [0, A] such that P[A > 0] = 1 and, for every 
M, N £ S d , u G n and y e B u 

(Here ^ are the usual Pucci extremal operators; see the next section). 
(F3) Uniform continuity and boundedness: for each R > 0, 

{F(-,-,uj) : oj G fi} is uniformly equicontinuous on Br x R d 

and 

esssup|F(0, 0, oj) | < +oo. 

(F4) Bounded moment of the ellipticity: the random variable A satisfies 

E [\- d ] < +oo. 

The main result. We now present the homogenization result, which for simplicity 
we state in terms of the Dirichlet problem 

(1.7) (f( D *u^)=0 inU, 

I u £ — g on dU. 

Here U C R d is a bounded Lipschitz domain and 3 G C(dU), and the PDE is to 
be understood in the viscosity sense (c.f. [10|,|6[). By modifying our argument in a 
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very minor way (only small changes in part three of Section |4]), we may homogenize 
any other well-posed problem involving F, including parabolic equations like 



ut + F^u, -,ujj = 

with appropriate boundary/initial conditions. 

Note that by (Fl) and (F2), for each e > 0, the equation (|1.7[) is uniformly 
elliptic with probability one. Indeed, if we take {B(xj, 1) : 1 < j < k} to be a 
finite covering of e~ x U , then its ellipticity is bounded by the random variable 

A sup A _1 (t , / e w) , 

l<j<k 

which is almost surely finite by (F2). See (|2.1I) below. As a consequence, (|1.7I) is 
well-posed and has a unique viscosity solution u e = u £ (x,uj) belonging to C(U). 

The main homogenization result is: 

Theorem 1. Assume (Fl), (F2), (F3) and (F4). Then there exists an event 
f2i G U of full probability, a positive constant < Ao < A which depends only 
on d, A and E[A~ d ] and a function F : E> d — > R which satisfies 

n AM — N) < F(M) - F(N) < 7t Q A M ~ N ) 

such that, for every w € fli, every bounded Lipschitz domain U C R d and each 
g G C(dU), the unique solution u £ — u £ {x,uj) of the boundary value problem (|1.7[) 
satisfies 

lim sup \u e (x,Lo) — u(x)\ = 0, 

where u G C(U) is the unique solution of the Dirichlet problem 

= inU, 
on dU. 



(1.8) 




A brief literature review. The homogenization of linear uniformly elliptic equa- 
tions in random media originated in the work of Papanicolaou and Varadhan (25|.|26| 



and Kozlov [23|,l24|]. Later, Dal Maso and Modica [ill. 1 1 21] obtained stochastic ho- 
mogenization results for nonlinear equations in divergence form and convex varia- 
tional problems. The homogenization of uniformly elliptic, nonlinear equations in 



nondivergence form was first considered in the periodic setting by Evans 15[ and 
much later by Caffarelli, Souganidis and Wang [8j in random media. In contrast to 
the divergence form case (c.f. [9(), little seems to be known about the homogenized 
coefficients for nondivergence form equations, even in the periodic case, other than 
what is inherited from the uniform properties of the medium. As far as quanti- 
tative homogenization results, we mention the work of Yurinskii [291 ] and Gloria 
and Otto 17, [l8| for linear equations and Caffarelli and Souganidis 7} for fully 



nonlinear equations. 

The modern regularity theory for elliptic equations in nondivergence form began 
in the 1980s with the groundbreaking work of Krylov and Safonov, Evans, Caffarelli 
and others, and we refer to [6] and the references there for more. Recently, there 
has been some progress in obtaining Harnack inequalities and Holder regularity 
for certain nonlinear degenerate equations (see for example 0,0, HH). In these 
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works, the degeneracy of the equation is typically compensated in some way by 
dependence on the gradient. A typical model equation considered is 

(1.9) \Du\' ( F(D 2 u) = 0, 

where 7 > and F is uniformly elliptic. A solution u of (II. 9[) may only be irregular 
if \Du\ is small, and this allows to compensate for the degeneracy. This is a very 
different situation from the "naked" degeneracy of the equations considered here. 

Outline of the paper. In the next section we give some preliminary results and 
notation needed later in the paper and make some comments about our assump- 
tions. In Section [3] we develop the deterministic regularity theory. The proof 
of Theorem [1] is then given in Section |4] Finally, in Section [5] we construct an 
explicit example to show that the moment condition (F4) is sharp. 

2. Preliminaries 

In this section we present some background results needed in the rest of the pa- 
per, including the statements of the ergodic theorems we cite, some remarks about 
our model and some general remarks concerning viscosity solutions and semiconcave 
functions. We begin by reviewing the notation. 

Notation. The symbols C and c denote positive constants which may vary at each 
occurence and which typically depend on known quantities. We work in Euclidean 
space R d for d > 1. We denote the set of natural numbers by N := {0,1,...} 
and Q is the set of rational numbers. The family of bounded Lipschitz subsets 
of M. d is denoted by L. The open ball centered at y G M. d with radius r > 
is B r (y) := {x G M. d : \x - y\ < r} and we write B r := B r (0). If E C R d is 
a bounded Borel set, then E is its closure and \E\ is the Lebesgue measure of 
E. If / G L l (E), then the average of / in E is f E f(x) dx := l^]- 1 f E f(x) dx. If 
/ : E — > R, then we denote osce / := sup B /— mfg /. The characteristic function of 
a Borel set E is xe- We work with a probability space (0, 2r, P) as described in the 
previous section. The indicator random variable of an event A e Jis written 1^. 
We say that A G 3 is of full probability if F[A] = 1. The space of symmetric d-by-d 
matrices is E> d . If M, JV e S d , we write M > N if the eigenvalues of M — N are 
nonnegative. If x,y 6 M. d , then x®y denotes the d-by-d matrix with entries (xtyj). 
The trace of M e § d is tr(Af). Recall that any M £ E> d can be uniquely expressed 
as a difference M = M + - M_ where M+M_ = and M+, Af_ > 0. In particular, 
if r 6 R, then we write r + := max{0, r} and r_ := (— r)+. The Pucci extremal 
operators are defined for < /1 < A and M G § d by 

T+ A (M) := -/xtr(M+) + Atr(M_) and V'^M) := -Atr(M+) +/xtr(M_). 

The elementary properties of the Pucci operators can be found in [6|. Here we 
remark only that they are uniformly elliptic, CP^ A is convex and A is concave. 

Brief remarks concerning the assumptions. Note that the restriction of (F2) 
to y G B\ is merely for convenience, it may be extended to all y G R d by stationarity. 
Indeed, the combination of (Fl) and (F2) yields, for all y, z G M. d with \y — z\ < 1, 

(2- 1 ) y A(r,o,),A( M - N ) < F (M, y, u>) ~ F(N, y, u) < 3^ (t , w)iA (M - N). 
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In light of (|2.1j) . it is convenient to abuse notation by writing X(z,u>) — X(t z lo). 
Note also that due to (F3) we may suppose that 

(2.2) {X(-,u}) : weSl} is uniformly equicontinuous on Mr. 

Otherwise we simply redefine A to be the largest quantity which satisfies (F2), 
which then satisfies (|2.2[) by (F3). The operators on the leftmost and rightmost 
side of (|2.ip are the minimal and maximal operators, respectively, which satisfy con- 
ditions (F1)-(F3). In particular, since X(-,uj) > 0, our equation is locally uniformly 
elliptic in the sense that, almost surely, imV A(-,o;) > for each V € £■ 

Using ergodicity, we may improve the second part of (F3) to 
sup ess sup \F(0, y,ui)\ < +oo. 

Using then the continuity of F and intersecting the event on which the latter holds 
over all the rational points of R d , we obtain 

ess sup sup \F(0,y,u))\ < +oo. 

Applying also (F2), this yields, for Cq := esssup weS1 \F(0,y,uj)\ and all M G E> d , 

(2.3) ess sup sup \F(M,y,u)\ < C + Atr(M+ + M_) < C(l + |M|). 

wen yt£m d 

A brief remark concerning viscosity solutions. All differential inequalities in 
this paper are to be interpreted in the viscosity sense (c.f. [(3, [HI). We remark that, 
while it is not obvious- in fact, it is equivalent to the comparison principle- we have 
transitivity of inequalities in the viscosity sense (see 0, Lemma 3.2]). For example, 
if V € L and u, v € C(V) satisfy 

F(D 2 u,y,uj) > and F(D 2 u,y,uj) < in V 

then formally it follows that for w := u — v we have 

(2.4) < F(D 2 u 7 y,u J )-F(D 2 u,y 7 u J ) < T+ {x u) A (D 2 w) . 

We emphasize that we may also deduce fj, ^ A (D 2 w) > in the viscosity sense, 
and make other similar formal deductions rigorous, using 0, Lemma 3.2]. 

Pointwise notions of twice differentiability and G > 1 . We require the follow- 
ing pointwise regularity notions. We say that u £ C'(B(0, 1)) is twice differ entiable 
at x S Bi if there exist (A, p) G S d x R d such that 

limsup sup r~ 2 \u(y) — u(x) — p ■ (y — x) — ^(y — x) ■ X(y — x) | = 0, 

r->-0 y eB r {x) 

in which case we write D 2 u(x) := X and Du{x) :— p. We also say that u is C 1,x 
on a set E C B(0, 1) if u is differentiable at each point of E and 

\u(y) — u(x) — Du(x) • (y — x)\ 
sup sup — — .1 ' yy < +oo. 

xeE V eB 1 \x-y\ 2 

A function u is semiconcave if there exist k > such that the map x i— > u(x) — k\x\ 2 
is concave. In this paper we rely many times on the observation that, for any a > 0, 
a semiconcave function is C ' on the set of points at which it can be touched from 
below by a C 2 functions with Hessian bounded by a. Moreover, by Rademacher's 
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theorem and the Lebesgue differentiation theorem, any C 1 ' 1 function on a set E is 
twice differentiable at (Lebesgue) almost every point of E. 

Infimal convolution. We recall a standard tool (c.f. [HI 0] for details) in the 
theory of viscosity solutions. If u G C(i?i) then its infimal convolution is denoted 

by 

u £ (x) := inf [ u(y) + -\x - y\ 2 J . 

The function u £ is more regular than u and, in particular, is semiconcave. It is a 
good approximation to u in the sense that u e —> u locally uniformly in B± as e — > 0. 
Moreover, if /, A G C(-Bi), A > and 

^(x),a(^)>/ m Si, 

then there exist sequences of functions X' s , f' e G C(B\) which converge locally uni- 
formly to A and /, respectively, as e —> 0, such that u e satisfies 

^KAPV) > f'e in B l-r c 

where r e — > as e — > 0. We refer to for details. For us, the principle utility of 

these approximations is the semiconcavity of u e . If u e can be touched from below 
by a smooth function <p at some point z £ B\, then u £ is C ' at z, with norm 
depending only e and \D 2 ip(z)\. See 0, Theorem 5.1]. 



Statements of the ergodic theorems. We next recall the two versions of the 
(multiparameter) ergodic theorem used in this paper. The first is nearly a conse- 
quence of the second, but since it is simpler we give it separately. A nice proof can 
be found in Becker [4| . 

We emphasize that the assumptions on the probability space (fi, 5F, P), in par- 
ticular (|1.5|) and (|1.6j) . are in force. Recall that L denotes the set of all bounded 
Lipschitz subsets of M. d . 

Proposition 2.1 (Wiener's ergodic theorem). Let f S L 1 (J7). Then there exists a 
subset 51o £ J of full probability such that, for every lo S fio md V G L, 

(2.5) lim / f(r v uj) dy = E[f] . 

In particular, the map y i— > f(T y u)) belongs to Lj l oc (M ) for every w G £Iq. 

The version of Proposition ^. ll proved in [J] actually requires V to be star-shaped 
with respect to the origin. As is well-known, this restriction may be removed as 
follows. First we notice that the conclusion holds for any cube V = Q with sides 
parallel to the coordinate axes, since any such cube either contains the origin or 
has the property that, for some larger cube Q, both Q\Q and Q are star-shaped 
with respect to the origin. Since it holds for such cubes, it holds for an arbitrary 
finite disjoint union of them and hence any V G £ by approximation. 

We next state the multiparameter subadditive ergodic theorem of Akcoglu and 
Krengel [l| as modified by Dal Maso and Modica 



I'M , which requires some further 



REGULARITY AND STOCHASTIC HOMOGENIZATION 



9 



notation. We denote by Uq the family of bounded subsets of M. d . A function 
/ : Ho — > R is subadditive if, 

fe 

f(A)<J2f(A J ) 

3 = 1 

whenever k e N and A,Ai,...,Ak € Uq are such that U* =1 Aj C A, the sets 
j4i, . . . , A\. are pairwise disjoint, and \ A \ Vfc_ x Aj \ — 0. Let M be the collection of 
subadditive functions / : Ho — > K which satisfy 

< f(A) < \A\ for every A E U . 

A subadditive process is a function / : fi — > M. It is sometimes convenient to write 
/(A, = /(w)(j4), in which case we have f(A, T y ui) = f(y + A, u>). 

Proposition 2.2 (Subadditive ergodic theorem). Let f : — > M be a subadditive 
process. Then there exists an event fio 6 2r of full probability and a constant 
< a < 1 such that, for every uj € fio md V G L , 

(2-6) Inn ^ = a. 

^ ' t^oo \tV\ 

This version of the subadditive ergodic theorem is [l], Proposition 1] in the special 
case that L is replaced by the family of all cubes, and we recover the general case 
by an easy approximation argument. 

3. Regularity in the macroscopic limit 

The classical regularity theory (as developed for example in Q) does not directly 
help us to homogenize (|1.4[) because, as e becomes small, the ergodic theorem 
guarantees that the set where (|1.4|) has very high ellipticity becomes dense. What 
we need are estimates which do not degenerate as e — s- 0, and for this it is necessary 
to revisit the regularity from the beginning. 

What the ergodic theorem ensures is that, almost surely in uj, for every > 0, 

(3.1) lim/ A- d (-,w) dx = E[\- d t {x< ^ } ] . 

Jvn{\<ti} \e ' 

In this section we develop a deterministic regularity theory for solutions of (|1.4j) 
which will be robust in the almost sure macroscopic limit e — > by virtue of (|3.ip . 

Since the random environment plays no role here, we drop dependence on uj. 
Throughout this section we consider a continuous function A : R d — > (0, A], and we 
study the regularity of subsolutions and/or supersolutions of the extremal operators 
^\(x) A m bounded Lipschitz domains V E £>. Our estimates must depend only on 
d, A and, for /i > 0, the quantities 

\- d (x) dx. 

'vn{\<^} 

As it is purely deterministic, the regularity developed here is of independent interest. 

The primary goal is to obtain an improvement of oscillation result on unit scales, 
giving us a modulus of continuity (and a Holder estimate in the macroscopic limit 
for solutions of (|l-4jl ). Our arguments are loosely based on the arguments in the 
classical regularity theory with some nice modifications due to Savin ^gfy, but 
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require several new ideas to overcome the degeneracy of the equation. For instance, 
"two important tools" are introduced in @, Section 4.1] which are used repeatedly 
in what has become the standard proof of the Harnack inequality. In our situation, 
neither of these tools can be applied in a straightforward manner. 

First, showing that an appropriate barrier (or "bump") function exists- which 
is easy in the uniformly elliptic situation (see [a, Lemma 4.1])- is a very nontrivial 
matter. Here we construct a barrier by touching a candidate function from below by 
translations of the singular function \x\~ a with a>l and then adapting the proof 
of the ABP inequality to show that the corresponding contact set would be too 
large if the function failed to be a barrier. Second, the measure-theoretic argument 
involving the Calderon-Zygmund cube decomposition must be altered due to the 
presence of "bad" cubes of high ellipticity, and we actually use an alternative new 
idea based on the Besicovitch covering theorem. 

The development is essentially self-contained and depends also on some novel 
uses of the area formula for Lipschitz functions, similar to the proof of the ABP 
inequality, and partially inspired by Savin [27l . 28 1. The main task of this section 



is to prove the following improvement of oscillation result. 

Proposition 3.1 (Decay of oscillation). There exists S > 0, depending only on d 
and A, such that if < /i < \ and 



f 

■IB 



\- d (x) dx < S, 



then there exist constants < r < 1 depending only on d, A, and \x, such that for 
all a > and u G C(B±) satisfying 

y l(x)A D2u ) ^ - Q and y A(x),A( I)2 ' lt ) < o *" B x , 

we have 

osc u < toscw + a. 

Bx/4 Bi 

Proceeding with the proof of Proposition 13. 1[ we begin with three applications 
of area formula for Lipschitz functions (c.f. [l6|), which asserts that 



\f{E)\ = f \det Df(x)\ dx, 

J E 



for all Lebesgue measurable sets E C WL d and injective Lipschitz maps / : E — >• M. d . 

The first is the Alexandroff-Bakelman-Pucci (ABP) inequality. The version we 
give here is not new: it is actually a corollary to the proof of @, Theorem 3.2]. 
We include a proof below both for completeness and in order to introduce the 
style of argument we use below, in a more complicated form, to obtain the barrier 
function. The argument here is much simpler than the one in 0, which is due to 
the observation that a semiconcave function is necessarily C 1,1 on the set where it 
can be touched from below by a plane. 

Proposition 3.2 (ABP inequality). Ifu,f £ C(i?2) satisfy 
1u>0 onB 2 \B u 
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then 



u_(0) < ( — / A "|.r )/;(, )</,- 

w/iere 



l-Bll J{T u =u} 



T u (x) := sup inf (p ■ (x - y) - U-(y)) 
is convex envelope of —U- := min{0,it}. 

Proof. By approximating u by its infimal convolution, we may assume that u is a 
semiconcave (it is straightforward to see that the limsup of the contact sets for u e are 
contained in the contact set for u). Also assume that a := u_(0) > 0, and observe 
that, for every p e B a , there exists z(p) S -B2 such that u(z(p)) < and the map 
x 1 ^ — m_(.t) —p-x attains its infimum in Bi at a; = Note that we can arrange 
for z : B a —> B2 to be Lebesgue measurable by choosing, say, lexicographically 
among the minimizers closest to the origin. Since u is semiconcave and can be 
touched from below by a plane on A :— z(B a ), it is C ' on A. In particular, ~z has 
a Lipschitz inverse p(z) :— Du(z). Moreover, u is twice differentiable at Lebesgue 
almost every point of A by Rademacher's theorem (c.f. [ijj]) and the Lebesgue 
differentiation theorem. At such regular points z £ A, we have that D 2 u{z) > 0, 
since u can be touched from below by a plane at z, and then the supersolution 
inequality gives 

"/(*) < y Mz),A (DM-A) = -Hz) tr {D 2 u{z)) ■ 
We deduce that, for almost every point z of A, 
(3.2) < D 2 u(z) < \~ 1 (z)f + (z)I. 

An application of the area formula yields 

|B 1 |a d =|B„|= / |det-D^)|dx < f \- d (x)ff{x)dx. 

J A J A 

Since A C {T u = u}, we obtain the proposition. □ 

Using a more sophisticated version of the above argument, we next construct 
the barrier function, which below plays a critical role in the proof of Lemma 13.61 
below, similar to that of the "bump" function in 6|. It is also needed in the next 
section in proof of Theorem Q] to verify that the limit function satisfies the Dirichlct 
boundary condition. 

Lemma 3.3. For each < r < ^, there exists S > 0, depending only on d and A, 
such that if < fx < 5 and 

f \- d (x)dx<5r d , 
J B!n{A<^} 

then there exists a constant f3 > 0, depending only on d, A, r, and fj,, such that for 
each u £ C(B\ \B r ) satisfying 

< u > on dBi, 

u > /? on dB r , 

we have u > on B\_ r \ B r . 



12 



S. N. ARMSTRONG AND C. K. SMART 



Proof. By approximating u by its infimal convolution, we may assume that u is 
semiconcave (alternatively, by the comparison principle we may replace u by the 
solution of the Dirichlet problem, which is C 2 by the Evans-Krylov theorem since 
Sa'w.a is convex). 
Define 

2(d-l)A , n /4 

f and $ — ' 



The main observation is that if u is negative somewhere in B\_ r \B r , then it can be 
touched from below by small translations of the singular function 4>{x) := 2 a \x\~ a . 
Indeed, suppose that u(xq) < for some xq £ B±- r \ B r . Then, for every y € B r / 2 , 

(3.3) inf (u(x) - 4>{x - y)) > -2<* 1 - T - ° > u(x ) - <f>{x - y) 
and, by our choice of /3, 

(3.4) u(x) > (3 > (f>(x — y) for every x £ dB r . 

Hence for every y G -B r /2, the map x i— > u(x) — <f>(x — y) attains its infimum at some 
point ~z{y) S B\\B r . We can easily arrange for the function z : B r / 2 — > B\ \ B r to 
be Lebesgue measurable. To obtain a contradiction, we eventually apply the area 
formula to the inverse of z. Most of the rest of the argument is concerned with 
showing that z has an inverse y and estimating | det Dy\. 
The differential inequality for u at z = z(y) yields 

- 1 < y + x (z) ^ D2 ^ z y)) - a ( a + l ) T \ z - y\~ [a+2) (f - x{ ~ z )) ■ 

Using that 2 a \z - y\-( a + 2 ) > i and a(a + 1) > 8/n, we obtain 

(3.5) \(z) < fi + - T ±— < [i. 
We conclude that 

(3.6) A := z(B r/2 ) C B 1 n {A < /i}. 

Since it is semiconcave and |D 2 </5| is bounded in M. d \ B r / 2 , we see that u is C 1,1 
on A. In particular, it is differentiable at each point of A and, by Rademacher's 
theorem and the Lebesgue differentiation theorem, twice differentiable at Lebesgue 
almost every point of A. For each y G B r / 2 , 

Du(z(y)) = D4>(z{y) -y) = -a2 a \z(y) - y\- (a+2 \z(y) - y). 

From this we see that ~z has Lipschitz inverse y : A — > B r / 2 given by 

y(z) :=z+ {a2 a )- 1 \z - y\ a+2 Du{z) = (a2 a ) 1/(a+1) \Du{z)\- {a+2)/{a+1) Du(z). 

At points z G A where u is twice differentiable, we have 

(3.7) - a2 a \z - y{z)\-^ a+2 h < D 2 (f>{z - y{z)) < D 2 u(z) 
and, using the differential inequality, 

_1 ^ y \(z) A ( D2u ( z )) ^ Atr ( £,2m (^))_ - H x ) tr(D 2 u(z)) + . 
A rearrangement of the later gives, in light of (|3.7|) . 

ti-(D 2 u(z))+ < (dka2 a \z-y{z)\- {a+2 ^ + l) \- l {z) 
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and from this we deduce that 

(3.8) D 2 u(z) < CX{z)- 1 (a2 a \z - y(z)\- (a+2) + 1)1, 

where here and in the rest of the proof C > depends only on d and A. Combin- 
ing (13. 7[) and (13.81) . we obtain, at each point zeiat which u is twice differentiable, 

(3.9) \D 2 u{z)\ < C\-\z) (a2 Q |z - y{z)\- (a+ ^ + l) . 

Differentiating the formula for y and using the expression for Du(z(y)), we get 
Dy(z) - I 

= (a2 a )- 1 \z - y(z)\ a+2 D 2 u(z) + (a2 a )- 1 (a + 2)\z- y(z)\ a (z - y) ® Du{z) 
= (a2 a )- 1 \z - y{z)\ a+2 D 2 u{z) + (a + 2)\z - y(z)\- 2 (z - y(z)) ® (z - y(z)). 

Using (l3Jt and that 2~ Q |z - y(z)\ a+2 < 4, we obtain 

\Dy{z)-I\<C\-\z) + {a + 2). 

Noticing that Q3.6P implies that a + 2 < CX^ 1 (z) for every z € A, we at last deduce 

(3.10) \Dy{z)-I\<C\-\z) 
for Lebesgue almost every point of A. 

We may finally apply the area formula, using (13.61) , (13.10[) and the hypothesis of 
the lemma to conclude that 

(2- d \B 1 \)r d = \B r/2 \= f \det Dy(x)\dx 

J A 

< C [ \- d (x)dx < C [ \- d (x)dx < C5r d . 

J A J B 1 n{\<n} 

We get a contradiction if 5 > is sufficiently small, depending on d and A. □ 

Remark 3.4. By an easy modification of the above proof, we obtain a slightly 
improved estimate for the barrier. Given h > we can modify the choice of 5, j3 > 
and also select an r' G (r, 1) depending only on d, A, e, /Lt, such that the supersolution 
u satisfies u > (3 — h on B r > \ B r . We use this observation to verify the Dirichlet 
boundary condition for the limit function in the proof of homogenization below. 

The next lemma, which is inspired by (27l . Lemma 2.1], follows from another 
application of the area formula and a (much easier) variation of the above argument. 
It asserts that, if a supersolution can be touched from below by sufficiently many 
translations of a fixed parabola, then the L d norm of A -1 on the set of points at 
which the touching occurs cannot be too small. 

Lemma 3.5. Suppose that u G C(B\) satisfies 

?W^ 2u ) > -1 in*, 

and there is an open set V C Mr and a constant a > 1 such that the map z >— > 
u{z) + ^\z — y\ 2 attains its minimum in B\ for all y G V. Let W C B\ denote the 
set of points at which, for some y G V, this map attains its minimum. Then there 
exists a constant 5 > 0, depending only on d and A, such that 

[ \- d (x)dx > S\V\. 
Jw 
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Proof. By replacing u by u + a\x\ 2 and letting a — > 0, we may suppose that, for 
some small r\ > 0, W C Bi-,, and, for every 

(3.11) mm + ||z - y| 2 ) > inf (u(z) + %\z - y\ 2 ) . 

As in the previous two arguments, we may assume that u is semiconcave by inhmal 
convolution approximation. Indeed, due to p. lip , the set V is essentially unchanged 
by the infimal convolution approximation, while the set W is unchanged or even 
smaller. 

Select a Lebesgue-measurable function z : V — > B\ such that the map z >->• 
u(z) + %?\z — y\ 2 attains its infimum in B\ at z — ~z(y). The function u is C ' on 
A := z(V) and z has a Lipschitz inverse y given by 

y(z) :— z H — Du(z). 
a 

By Rademacher's theorem and the Lebesgue differentiation theorem, u is twice 
diffcrentiable at almost every point of A and —al < D 2 u at such points z, where 
the differential inequality also gives 

-1 < 7+ {z) K {D 2 u{z)) =Atr (D 2 u{z))_ - \{x) tr{D 2 u{z)) + . 
From this we obtain that, almost everywhere on A, 

(3.12) - al < D 2 u(z) <J + adAA(z)" 1 / < CA^ 1 (z) (1 + a) I. 
An application of the area formula for Lipschitz functions yields 

|V|= f \det Dy(x)\ dx ^Cil + a- 1 ) 11 [ ^ d (x)dx, 

J A J A 

from which we obtain the lemma, using that a > 1 and A C □ 

We now give the proof of the decay of oscillation. 
Proof of Proposition \3.1[ Now suppose that < /i < | and 

(3.13) / A- d (.x)da:< ( 5:=-i r |B 1/ 3|min{8-%,32-%}, 
JS2n{A< M } o-'M 

where S± is from Lemma 13.31 S% is from Lemma 13.51 and Nd is the constant from 
the Besicovitch covering theorem in dimension d. 

We first make a reduction by noticing that, to obtain the proposition for r := 
(1 — 1/fc), it suffices to consider v 6 C(-Bi) which satisfies 

(3-14) ^WOD 2 ^" 1 and P a"(*),aP 2 *) < 1 

and 

(3.15) oscv < 1 + osc v, 

Bi -B1/4 

and to show that osc^ v < k. Indeed, suppose u satisfies the hypotheses of the 
proposition and 

osc u > (1 — k ) osc u + a. 

Bi/4 Bi 

Then a < oscsj u, and so if we set v := fcu/ oscs 1 it, then we see that v satisfies 
(|3T4|) . ([335]l . and osc Bl v = k. 
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Define, for every n > 0, 

An ■= [x G .Bi : 3y G B u v(x) + ^\x - y\ 2 = Jrf + ||z - y| 2 ) J . 

In other words, A K is the set of points in B1/3 at which v can be touched from 
below by a paraboloid with Hessian — kI and vertex in B\. To prove the desired 
estimate on u, it is enough to show that 

(3-16) \A K KB X/S \ > §|B 1/3 | 

for some k > depending only on d, A and \i. Indeed, if we could show this, 
then using that — v satisfies the same hypotheses as v and applying (|3.16j) to both 
functions, we find a point x G S1/3 which can be touched from above and below by 
parabolas with opening n. That is, there exist x G -B1/3 and 2/1,2/2 G B\ such that, 
for all z E Bi, 

u(a;) + -|a;-yi| < u(z) < v{x) - -\x - y 2 \ + -\z - y 2 \ . 

This implies that \v(x) — v(z)\ < 8k/9 < k for all z G Bi, and thus oscsj^ d < 2k, 
which is the desired estimate. 

In order to obtain p.l6j) for some k > depending on the appropriate quantities, 
we observe first that (|3.15[) implies that A144 n i?i/3 ^ 0, since v can be touched 
from below in -B1/3 by the parabola — 144\x — y| 2 , where y G -B1/4 is such that 
v(y) = min^ v. We then repeatedly apply Lemma l3^6l below to obtain the desired 
result for k = 1440 n , where n := \\Bi\/r]~\ and 8, r\ > are given in the statement of 
the lemma. The proof of Proposition 13. II is now complete, pending the verification 
of Lemma 13.61 □ 



To complete the proof of Proposition 13.11 we must prove the following lemma. 
In the classical regularity theory, a result like this typically relies on the Calderon- 
Zygmund cube decomposition (as in the proof of (4.12) in [6]). Since we did not 
see how to adapt it, we instead use an alternative measure-theoretic tool, the Besi- 
covitch covering theorerrQ. The argument also relies on Lemmas 13.31 and 13.51 



Lemma 3.6. Let [i, v and A K be as in the proof of Provosition \ 3.1\ There exist 
constant 9 > 1 and i] > 0, depending only on d, A and /i, such that if k > 1, 
A K n B 1/3 + and \A K n B 1/3 \ < ||B 1/3 |, then \Ag K \ > \A K \ + 77. 

Proof. Consider the collection 23 of balls B r (x) C B 2 such that B r / 2 (x) C B\ \ A K 
and dB r j 2 (x) n A K ^ 0. Note that since A K n B x / S ^ and A K is closed, every 
point of -B1/3 \ A K is the center of some ball in 23. According to the Besicovitch 
covering theorem, we may select a countable subcollection {B rk (xk)}keN ^= 23 that 
covers B 1 / 3 \A and such that each point x G B 2 belongs to at most A^^ balls. 

We say that the ball B rk (xk) is good if 

\- d {x)dx < | J B rfc (x)|min{8-%,32- ti <52} 

B r (x k )n{\<ti} 



^Luis Silvestre has since pointed out to us that the Calderon-Zygmund decomposition argu- 
ment in |6|] may indeed be suitably modified to prove Lemma 13.61 and that, alternatively, the 
Vitali covering theorem may be used instead of the Besicovitch covering theorem with only small 
modifications to the argument we give below. 
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and set G :— {k G N : B Tk (xk) is good}, where 81,82 > are as in the proof of 
Proposition 13. II We claim that at least half of the measure of -B1/3 \ A K consists of 
points which belong to good balls, i.e., 



(3.17) 



U B rk (x k ) 



>hB 1/3 \A K \ > hB 1/3 \ 



Indeed, if (|3.17p were false, then J2kgG \Br k \ > 5 \^i/a \ Ac| > ^\B 1 / 3 \ and so 



J B 2 n{\<^} ly d k g G JE 



\- d (x)dx 



' B r (x k )n{\<n} 



> ^-min{8- d 8 u 32-%} £ \B r J > -L m m\B 1/3 \ 32-%} , 

d k$G d 

which contradicts (|3.13p . Therefore, in light of the Besicovitch covering, it is enough 
show that \B,. k (x k ) n Ag K \ > rj\B r . k {xk)\ for each good ball B,. k (xk) and some 
constants 9, r\ > depending only on d, A, and /x. 

Fix a good ball B r (x) := B Tk (xk) and choose z\ G dB r / 2 (x) R A K . By the 
definition of A K , there exists y\ G B\ such that 



(3.18) 



w 0i) + 2 ki - yi I 2 



inf (v(z) + -| 



z-yi\ 



Applying (a properly scaled) Lemma [3.3[ since B r (x) is good there exists (3 > 1, 
depending only on d, A and [A, such that the solution w of the Dirichlet problem 



< 



( K (x) A d2w ) 

w = 

t iu = fir 2 



-1 in 5 r (x) \ B r/S (x), 
on dB r (x), 
on dB r/s (x), 



satisfies 10 > in B r / 2 {x) \ B r / S (x). Clearly w < fir 2 in B r \ B r / S (x) by the 
maximum principle. Observe that the function 



satisfies 
(3.19) 



tp{z) := (dAn + 2)w - -\z - yi \ 2 , 



^.aPV) < -2 in B r (:r)\B r/8 (x). 



The comparison principle implies that the map z t— > w(z) — V 5 ( z ) attains its 
infimum in B r (x) \ B r / 8 (x) at some point z — z 2 G dB r (x) U dB r /g(x). However, 
Z2 G dB r (x) is impossible due to (|3.18j) and the fact that w = on dB r (x) and 
ui(zi) > 0. Hence z 2 G dB r / s (x), and so <£(-z 2 ) = — f \zz — yi\ 2 + (dAK + 2)fir 2 . We 
therefore obtain, using that w > in B r / 2 (x) \ B r /g(x), 



(v(z) + ^\z- yi \ 2 ) > v(z 2 ) + I \z 2 - y x | 2 - (dA + 2/k) nfir 2 . 
Using (13.18[) . z\ G B r / 2 (x) and k > 1, we obtain 



mm 

z£B r / 2 \B T / s 



(3.20) inf (v(z) + ^\z- yi \ 2 ) >v{z 2 ) + ^\z 2 - yi \ 2 -(dA + 2)Kfir 2 . 

zGB 1 \ Z J Z 
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It follows that, if we set 7 := 16/3(dA + 2) + 1, then for every y 2 G B r / S (x), the 
function 

t/>(z) := + -|z - yi| + — |z - 2/2 

satisfies ^(22) < rninsi\B r / 2 (a:) "0 an d therefore must attain its infimum in i?i in 
B r / 2 {x). Consider the function z : B r /${x) — > B\ given by ~z(y) = (j/i +7J/)/(l +7). 
Observe that for some constant a G R, 

||z - ^l 2 + ?£\z V 2? = ^^\z - 2(2/2) | 2 + a for all z G R d . 

It follows that the map z 1— > v(z) + -rr\z — 2(7/2) | 2 attains its infimum in B\ at some 
point of B r / 2 (x). Since 7 > 1 and thus 7/(7 + 1) > i, we deduce that 

z(B r/8 ) > 2- d \B r/a (x)\. 

Applying Lemma 13.51 and using that B r (x) is a good ball, we obtain 

2- d S 2 \B r/8 (x)\ < [ \- d {x)dx 

J B r/2 (x)nA ( ^ + 1)lt 

< 32-%\B r (x)\ + ^ d \B r/2 {x) n A (7+1)K |, 
which implies \B r / 2 {x) Ci A( 7+1 ) K | > fi d (l — 2~ d )5 2 \B r / s (x)\, as desired. □ 

4. HOMOGENIZATION 

The proof of homogenization follows the approach of [1] , although we have re- 
organized the argument for clarity and simplicity as well as to accommodate the 
modifications required to handle the non-uniformly elliptic case. The reader inter- 
ested only in the uniformly elliptic case can find here a concise proof of the result 
of . The strategy relies on an application of the subadditive ergodic theorem to 
a certain quantity involving the obstacle problem. The proof has three steps: 

(1) Identifying F: by applying the subadditive ergodic theorem to the Lebesgue 
measure of the contact set of a certain obstacle problem, we build the 
effective operator F. 

(2) Building approximate correctors: with the help of the effective regularity 
results, we compare the solutions of the obstacle problem to the solution 
of the Dirichlet problem with zero boundary conditions and show that the 
latter act as approximate correctors. 

(3) Proving convergence: using the approximate correctors, the classical per- 
turbed test function method allows us to conclude. 

Step one: identifying F via the obstacle problem. Following Q, we intro- 
duce, for each bounded Lipschitz domain V G £, the obstacle problem (with the 
zero function as the obstacle): 

f mm / F(n 2 " 

(4.1) 



mm {F{D 2 w,y,uj),w} = in V, 
w = on dV, 



Some important properties of (|4.1[) are reviewed in Appendix [A] It is well-known 
that (|4.1[) has a unique viscosity solution, which we denote by w = w(y,cu ;V, F). 
We often write w = w(y,ui]V) or simply w = w{y,uj) if we do not wish to display 
the dependence on F or V . 
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The set Q(V,to) := {y £ V : w(y,uj; V) — 0} of points where w touches the 
obstacle is called the contact set. We write C(V, ui ; F) is we wish to display the 
dependence on F. The Lebesgue measure of this set is an important quantity, and 
we denote it by 

(4.2) m(V,w):=|e(V,w)|. 

We check that m satisfies the hypotheses of the subadditive ergodic theorem (Propo- 
sition I2.2[) . First we observe from the monotonicity of the obstacle problem, 
see (IA.4I) . that for all V, W £ £ and wgO, 

(4.3) VCW implies that G(W, uj) n V C Q(V, uj). 

Immediate from (|4.3I) is the subadditivity of m. That is, for all V, Vx, ■ ■ ■ , Vu £ £ 
anr-Vi that * Jj =1 Vj C V, the sets Vi, . . . , 14 are pairwise disjoint and |V\L$ =1 V}| = 0, 



we have 

k 

(4.4) m{V,u>) <^2m(Vj,ui). 

According to (Fl), m is stationary, that is, 

m(V, r y a;) = m(y + V, uj) 

for every y GM. d and We may easily extend m to Uo by defining, for every 

A £ 1X , 

m(A, uj) := inf {m(V, w) : and AC.V}. 

This extension agrees with m on £ by (|4.3|) and it is easy to show that the subad- 
ditivity and stationarity properties are preserved. 
We now obtain: 

Lemma 4.1. There exists an event ^2 £ 3" of full probability and a deterministic 
constant m £ R smc/i i/ia£ ; /or every uj £ and Lipschitz domain V C R d , 

(4.5) lim -^m(iV,a;) = to|V|. 

Proof. In light of the remarks preceding the statement, the lemma follows from the 
subadditive ergodic theorem (Proposition 12. 2[) . □ 



For clarity we write m = m{F) to display the dependence of m in Lemma 14.11 
on the nonlinear operator F. 

We are now ready to define the effective nonlinearity: 
(4.6) F(0) := sup {a £ R : m(F - a) > 0} . 

We extend this definition to all matrices in the obvious way. For each N £ S d we 
denote Fjv by 

F N (M,y,u>) := F(M + N,y,ui) 
and then we set, for each M £ § d , 

:=F M (0). 

To check that F is well-defined and finite, we first observe that, by (|A.2[) and (IA.6I) . 
inf F{0, y, uj) > implies that 6(1/, w) = V 
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and 

sup F(0, y, to) < implies that G(V, lo) = 0. 
yev 

Using (F3) and the remarks in Section [5J it follows from these that 

(4.7) essinfF(Af,0,w) <F(M) < ess sup F(M, 0,w). 

The monotonicity of the obstacle problem implies that a t— > m(F~a) is a decreasing 
function and thus m(F — a) > for a < F(0) and m(F ~ a) — for a > F(0). 

It is immediate from the comparison principle for the obstacle problem that, if 
Fi and F 2 are two operators satisfying our hypotheses, then 

(4.8) sup ess sup {F x (M, 0, lo) - F 2 (M, 0, w)) < implies F x < F 2 . 
Aies d wesi 

It is even more obvious that adding constants commutes with the operation F i— >■ F. 
From these facts a number of properties of F are immediate, the ones inherited from 
uniform properties of F. A few of these are summarized in the following lemma. 

Lemma 4.2. For every M, N G S d such that M < N, we have 

(4.9) < F(M) - F(N) < A tr(JV - M). 

Moreover, if M H> F(M, 0,lu) is positively homogeneous of order one, odd, or linear, 
then F possesses the same property. 

Proof. Each of the properties are proved using the comments before the statement 
of the proposition. To prove (|4.9|) . we simply observe that, according to (Fl), for 
all (Y,y,cu) G § d X R d x Q, 

(4.10) F(M + Y,y,u) < F(N + Y, y, u) + A tr(iV - M) 



and then apply (|4.8|) . It is obvious that F inherits the properties of positive homo- 
geneity and oddness from F, and linearity follows from these. □ 

Observe that (|4.9p asserts that F is degenerate elliptic. If F were uniformly 
elliptic, i.e., A -1 G L°°(f!,), then it follows from an argument nearly identical to 
the one for (j4.9[) that F is uniformly elliptic. In our more general situation, i.e., 
A -1 G L d (fl), the operator F is uniformly elliptic as well, but the proof is more 
complicated. We postpone it until the next subsection, since it is more convenient 
to deduce it as a consequence of Proposition 14.41 

Wc next show that, in large domains, the contact set has nearly constant density. 

Lemma 4.3. For every lo G ^2 and V, W G L with W C V, 
(4.11) hm — — — — = m. 

V ' t-yoc \tW\ 

Proof. Let U := V \ W G £> and fix u> G Vl 2 . Observe that (|4.3p gives 
(4.12 l msup^ — < Hm — - — =m 

1 ' t^oo \tW\ "t^oo \tW\ 

and, by the same argument, 

limsup 7——. < to. 

t^oo tU 
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Therefore, 

(4.13) liminf Ig^f^l = liminf R tV \ = ^ ° ^ 

Combining (j4~12|) and (|4~T3|) yields P~TT]t . □ 

Step two: building approximate correctors. The next step in the proof of 
Theorem Q] is to show that, in the macroscopic limit, the obstacle problem controls 
the solution of the Dirichlet problem 

j F{D 2 v,y,uj) = mV, 
) v = on dV. 



(4.14) 



As before, V G £ is a bounded Lipschitz domain and we write u = uj ; V, i* 1 ). 
The following proposition is the focus of this subsection. 

Proposition 4.4. There exists an event il 3 G £F o/ /uZZ probability such that, for 
every uj G 3) M G S d and V G £, 

(4.15) lim 4 sup w ; tV, F M - F(M)) I = 0. 
t^oo ^ y£tv 

Before we give its proof, we remark that Proposition 14.41 is a special case of 
Theorem [TJ We can see this by fixing U G £>, defining 

v e {x,uj) :=e 2 v { ~,uj;~U,F -F(0) 
\e e 

and then checking that v E (■, uj) is the unique solution of the boundary- value problem 

(4.16) f(dV,^)=P(0) inU, 

v £ = on dU. 



The conclusion of Proposition 14.41 then asserts that 
(4.17) v £ — > uniformly in U as e — >• 0, 

which is consistent with Theorem [T] since the zero function v = is obviously the 
unique solution 

j F(D 2 v) = F(0) in U, 
1 v = on 9Z7. 



(4.18) 



As we show in the next subsection, Proposition 14.41 actually implies Theorem [T] 
This is because, for large R > 0, the function £(y) := v (y,cu ; Br, Fm — F(M)j is 
an "approximate corrector" in in the sense that it satisfies the equation 

(4.19) F(M + D 2 ^, y, uj) = F{M) in Br, 

and is "strictly subquadratic at infinity" (i.e., satisfies (|4.15l0 . This is precisely 
what is needed to implement the perturbed test function method. 
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Proof of Proposition \4-4\ According to the ergodic theorem, there exists an event 
O4 G SF of full probability such that, for every ui E ^4, £ £ and rational q G Q 
with 9 > 0, 

(4.20) lim/ \- d {y,uj)dy = E[A- d l 
and 

(4.21) lim/ A- d (y,a;)x { A< 9 }(y)dy = E[A- d l { A< 9 }]- 

Note that, according to the ABP inequality (Proposition 13.21 properly scaled), for 
every uj G CI4, 

(4.22) lim lim sup ^ sup —r\v(y,uj;tV,F)-v(y,uj;tV ) F + a)\ = 0. 

We now define := 0,2 D H4, where ^2 is given in the statement of Lemma T4. II 
We first show that, for all uj G 3 , V G L, and M G § d 

(4.23) liminf 4 nrf w ; tV, F M - F(M)) > 0. 

We may assume that M = by replacing F with F_m and that -F(O) = by 
replacing F by F — -F(O). By (|4.22[) . we may also suppose that rn(F) — by 
considering F — a for a > and then sending a — > 0. Set 

if := esssup(i ;l (0, 0,w)) + = ess sup sup (F(0, y,u;)) + . 
weo wen a gR d 

According to (|4. 1 5[> and (|A.5|) . for every t > 0, the function it :— w(-,uj ;tV, F) — 
v(-, uj ; tV, F) satisfies 

^x(y,u)A D2u ) ^ K Xe(tv,.) in tV 
and u — on ^(iU). Using that u> > 0, the ABP inequality (Proposition 13.21 
properly scaled) and (|4.5p . we obtain 

(4.24) lim sup — sup — v(y, uj tV, F) < lim sup — sup u{y) 

t-¥oa t y£tV t— S-OO t y£tV 

<CK\imsup(J- \~ d (y,uj)xe(tv,u 1 )(y)dy) 

t^oo \JtV J 

To estimate the integral on the right, we observe that, for each k G N, 

^ d {y^)Xe(tv^){y)dy < ( k d \m(tV, u)\ + / \- d (y, u)x{\<i/k}(y) dy ) . 

\ Jtv J 

Divide this by \tV\ and pass to the limit t — > 00 using (|4.21[) to obtain 

(4.25) limsup / \- d (y, w) X e(tv,u)(v) d V ^ k d m(F) + E[\- d t {x<1/k} ] . 

f->oc JtV 

Since m(F) = 0, we may send k — > 00 and combine the resulting expression 
with (|4~24) to obtain (|3T2"3]) . 

To complete the proof, we show that, for every <j G ^3, 

(4.26) limsup ^ sup v(y, uj ; iV, F M - F(M)) < 0. 
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As above, we may suppose that M = and F(0) = 0. We may also assume that 
Tn(F) > 0, by considering F + a for a > and then sending a — > 0, using (|4.22l) . 
Since v < w, it suffices for (|4.26[) to show that 

(4.27) limsup^ sup w(y,uj;tV, F ) = 0. 

Furthermore, by the monotonicity of the obstacle problem it suffices to show that 

(4.28) limsup^- sup w (y,uj ;tB 2 R, F) = 

t— >oo t y£tB R 

where R > 1 is large enough that V C Br. Fix r > and observe that, 
by Lemma 14.31 and an easy covering argument using m(F) > 0, there exists 
T > sufficiently large such that, for every t > T and x £ Br, the function 
w(-,u) \tB2R, F) vanishes in B(tx,tr). We therefore have, for every t > T and 
x £ B R , 

(4.29) hw(tx,u;tB 2R ;F)\ < osc - s w(;u;tB 2R ;F), 

t B(tx,tr) t z 

We prove (I4.28P by showing that the lim-sup of the right-side of (|4.29l) . as t — > 00, 
is o(l) as r — > 0. For this we rely on Proposition l3.il 

Notice that (IA.5I) . (|4.20l) and the ABP inequality (Proposition yield, for 
t > sufficiently large, the bound 

(4.30) i sup \w(-,u;tB 2 R,F)\ < CKR 2 , 

t tB 2R 

where C depends only on d, A and E [A~ d ] . Select /i > such that 

E [ A ~ dl U<M}] < ^ 

where S > 0, 8 £ Q is as in Proposition 13. II By (|4.21l) . and making T > larger, if 
necessary, we have that for all t > T, r < r' < R and x £ Br, 

^ d {y,u)x{\ <P ,} dy < S. 

lB tr ,(tx) 

To see this, consider a finite covering {B s (xi)} of Br by balls of radius s = 2~ k R 
for some k £ N. According to (|4.21[) . for sufficiently large t, the average of 
\~ d (y,ui)x{\<(j,} in each of the balls B2 S (xi) will be less than A~ d 8. But every 
ball B r > (x), with s/2 < r' < s and x £ Br, is contained in one of the balls B2 S (xi). 
Since 4r' > s, this yields 

(4.31) / X- d (y,u)x{x<»}dy<4 d 'f A^(y, W ) X{ a <ai} dy < S. 

JB tr ,(tx) J B 2s (xi) 

Repeating this covering argument for k — 0,1,2, ... , |~log 2 (i?/r)] and making T > 
larger, if necessary, we obtain (|4.31[) for every t>T,r<r'<R and x £ Br. 

Iterating Proposition 13. II using (|4.30[) as well as (|A.3|) and (|A.5|) . we obtain, for 
every x £ Br and t >T, 

(4.32) osc ^w(;Lu;tB 2R ;F) < Cr^ 

for some constants 7 > and C > which may depend on d, A, E [A~ d ] , /i, K 
and R, but do not depend on r. Combining this with (14.291) and sending t — > 00 
and then r — >• 0, we obtain (I4.27[) and thus the proposition. □ 
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We conclude the second part by showing that F is uniformly elliptic and giving 
an estimate of its ellipticity. The proof is based on Lemma 15751 and Proposition ^. 41 

Proposition 4.5. There exists Ao > depending only on d, A and E[A~ d ] such 
that, for all M,N G § d ; 

(4.33) 3>- n A (M -N)< F(M) - F(N) < T+ A (A/ - N). 

Proof. Select matrices M, N £ E> d such that M > N. Fix w £ ^ and dehne, for 
each e > 0, 

V e (x) := e 2 v (-, u ; -B lt F M - F(M)) -s 2 v(-,uj; -B U F N - F(N)) 



J V e s I 



+ -x ■ (M - N)x. 



It is easy to check that V satisfies the inequality 

(4.34) y+ {i u)A (D 2 V e )>F(M)-F(N) in B x . 

According to Proposition 14. 4\ 



(4.35) V e (x) ->■ -x ■ (M — N)x as s -> uniformly in B\, 

Suppose that M — N has a largest eigenvalue a > with corresponding normal- 
ized eigenvector £ G R d , |£| = 1 so that 

(4.36) a£ ® £ < M - N < al 

Fix (3 > and, for each y G M d , denote by z(y) G R d the (unique) point at which 
the map x i-> $(x, y) := ^x • (M — AT)a: + /3|a; — y\ 2 attains its (strict) global 
minimum on R d . The inequalities <f>(z(y),y) < $(0,y) and &(~z(y),y) < $(y, y) 
imply, using (|4.36[) . that 

\a{z{y) • C) 2 + 0\z(y) - y\ 2 < min ( \a\y\ 2 , /%f 



In particular, we notice that \y\ < 1/3 implies that \z(y)\ < 2/3 and 



(4.:-!7) T(,y). £ 



1/2 



Applying (|4.35[) , we deduce that, for sufficiently small e > and every y G B 1 / 3l 
the infimum in B\ of the map a; i-> 14 (x) + /?|x — y| 2 is attained in -B2/3 and any 
point z at which the minimum is attained satisfies 

(4.38) l-^K£ 1/2 

Let A := B 2/3 D {x : \x ■ £| < (/3/a) 1 / 2 } and observe that |^| < C(/3/a) 1/2 - If 
F{M) — F(N) > —2/3, then we may applying Lemma [3751 using (|4.20[) . to obtain 

c< |Si/ 3 | < Climsup f \(-,u) d dx = C\A\E[\- d ] < C ( ^ f E[\- d ] . 
6^0 Ja^ £ ' 1 \ a J 
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This is impossible if a > C(3E[\- d ] . Therefore, we may select Ao > small 
enough, depending on the appropriate quantities, so that by fixing (3 :— \ ad/2, 
the previous argument yields, for all M > N, 

F(M) - F(N) < -2/3 = -X ad = 9+ oA (aI) < T+ o jM - N). 

Recalling (I4.9[) . we also have, for every M > N, 

F(M) - F{N) > -Atr(N - M) = ^ . A (M - N). 

We have verified (TQ51) for all M,N G § d with M > N. To remove the latter 
restriction, fix M, N G § d , write 

F(M) - F(N) = F(M) - F(M - (N - M)_) 

+ F(M - (N - M)_) - F(M - (N - M)_ + (N - M)+) 
and observe by what we have shown above that 

F(M) - F(N) < y+ _ A ((JV - M)_) - _ A ((JV - M) + ) = ?+ iA (M - iV). 
Interchanging M and TV yields the first inequality of for all M,N €§ d . □ 

Step three: concluding by the perturbed test function method. By adapt- 
ing the classical perturbed test function method, first introduced in the context of 
periodic homogenization by Evans (TBj , we now complete the proof of Theorem [1] 
The test functions are perturbed by the approximate correctors constructed in 
Proposition 14 . 4IT he argument we present here is similar in spirit to the one given 
in Section 4 of [8], although a bit less complicated. 

Proof of Theorem [TJ Fix a bounded Lipschitz domain U £ £, g G C(dU) and 
an environment ljq G f^, where the event G 5" is given in the statement of 
Proposition 14.41 

We first argue that, for every x G U, 
(4.39) u{x) :— lim sup w e (a;, wo) < u{x). 

To show (|4.39p . we begin by checking that u{x) < g on dU. By approximation, 
we may assume that g = and that U is smooth (and in particular has the exterior 
ball condition). By dilation, we may also assume that F(0,-,w) < 1 and that 
U C B R/2 (0). Given y G 9LT, we may select B r (x) CR d \U such that B r (x) nU = 
{y}. Given h > 0, we apply Lemma 13.31 with the modification in Remark 13.41 
Using u>o G ^3, we may select f3 > and r' G (r,R — r) such that the solution 
tp e G C(B R \B r ) of 

< = /3 on 9B B 

[ ip £ = on as r , 

satisfies limsup e ^ < /i in Vn B r >(x). Since C/ C Br{x) and u 6 < on 9t/, the 
comparison principle implies that u e < tp s . It follows that 

limsup sup u £ (-,cjo) < h. 

e-i-0 yns,.,_ r (y) 

Since h > was arbitrary, we conclude that u < g on <9£/. 
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By the comparison principle, to prove (|4.39p it suffices to check that the function 
u{x) := limsupg^Q u e (x, ujq) satisfies, in the viscosity sense, 

(4.40) F(D 2 u) < in U. 

To verify (|4.40[) . we select a smooth test function <j) € C 2 (U) and a point xq G U 
such that 

x i y (w 0) (x) has a strict local maximum at x = xo ■ 

We must show that F(D 2 (j>(xo)) < 0. Set M := D 2 cj)(xo) and suppose on the 
contrary that 9 := F(M) > 0. 

Since the local maximum of u — 4> at xq is strict, there exists ro > such that 
B rQ (xo) C U and, for every < r < ro, 

(4.41) (u - 4>) (x ) > sup (u -(/)). 

dB,,(x ) 

We next introduce the perturbed test function 

<t> e {x) := <j>(x)+e 2 v(-,uj ;-B ra (x Q ),FM-F{M) 
\e e 

We claim that, in some neighborhood of Xq, (ff is a strict supersolution of the 
oscillatory equation at microscopic scale e. More precisely, we will argue that, for 
some suitably small < s < ro to be selected below (and which may depend on 0), 

(4.42) f(dV,|,w„) > \e inB s (x Q ). 

To check (|4.42[) we select a smooth test function ip € C 2 (B s (xo)) and a point 
ii e B s (xo) such that 

,t i — y (<jf — V) i x ) has a local minimum at x = x\. 

Using the definition of (ff and rescaling, we have 

v(y,u ; ~B ro (x ),F M -F(M)j - ^ (ip(ey) - <f>{ey)) 



has a local minimum at y = — . 

£ 



Using the equation for i>, we obtain 



F^M + fl^nl-fl^n),^^ - F(M) >0. 

Since </> € C 2 , we may make \M — D 2 <fi(xi)\ = \D 2 tfi(xo) — D 2 cj)(xi)\ as small as we 
like by taking s > small enough. Thus, in light of (F2), we may fix s > so that 

F(M + D 2 iP{x 1 )-D 2 <j ) {x 1 ),^,uo) -f(d 2 ^{ Xi ),^,u;^ \ < U. 

The previous two inset inequalities and 9 — F(M) yield 

(4.43) >\e. 

This completes the proof of (14.421) . 

An application of the comparison principle now yields 

(4.44) u £ (x ,ujo) - (tf{x Q ) < sup (u e {-, cj ) - £ ) = sup (u £ (-, u ) - (ff) ■ 

B s (x ) dB 3 (x a ) 
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Taking the limsup of both sides of ()4.44|) as e — > and applying Proposition 14. 4[ 
we obtain 

u(xq) — 4>(xq) < sup (u — 4>) . 

dB s (x ) 

This contradicts (|4.41[) and completes the proof that F(M) < 0, and hence of (|4.40[) . 
and hence of ()4.39|) . 

It remains to show that, for every x G U , 

liminf u £ (x,lj ) > u(x). 

E-»0 

This is obtained by mimicking the argument above with only very obvious modifi- 
cations, and so we omit the proof. □ 



5. Breakdown of homogenization and regularity for p < d 

In this section we show that the condition that the dih moment of A -1 is finite is 
sharp for both the homogenization and regularity results. It was previously shown 
in [l9( that homogenization may fail unless some moment condition is imposed. 
Here we construct a stationary-ergodic random environment (Q, J, P, r) and an 
operator satisfying (Fl), (F2) and (F3) as well as 

(5.1) E[A~ P ] < +oo for all p < d, 

but for which homogenization fails. The example is actually for the linear operator 
F(M,y,uj) := — A(y, lj) tr(M), and we demonstrate the failure of homogenization 
(and an "effective" ABP inequality and hence effective regularity) by showing the 
solution u £ of the Dirichlet problem 

(5.2) |-A(^)A^ = 1 inQ, 

I u e — on dQ, 

satisfies lim £ _>o u e (0, w) = +00 almost surely, where Q := [—1,1] . The uniformly 
continuous function A : M. d — >• (0, 1] we construct below satisfies a finite range 
dependence, so even this strongest possible ergodicity/mixing assumption cannot 
save homogenization without the moment condition (F4). The argument that it E (0) 
blows up as e — >• is a combination of the weak Harnack inequality and ideas from 
percolation. 

The idea is to build "traps" where, from the probabilistic perspective, the cor- 
responding diffusion will become stuck for long periods of time. We choose, for 
each positive integer fc, a random arrangement Pfc(w) C M. d of points so that 
^(y,u) < Afc :— l/(fclog(2 + fc)) in each ball of radius 1 with center in P/c(w). 
To be more precise, we select, for each k € N, a continuous function 8k(y) on R d 
which is equal to A& in B\, takes the value 1 in R d \Z?2 and satisfies Afc < 9k(y) < 1- 
Then we set 

A(y,w) := inf inf 9 k (y-x). 

The point configurations Pk are independent Poisson point processes (c.f. [Hj]), 
with intensity such that the expected number of elements of Pk n V is equal to 
a\V\k~ 1 ~~ d , where a > is parameter we choose below. Since the series YlkLi k~ 1 ~ d 
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converges, it follows that the number of points of U^_ 1 Pk is almost surely locally 
finite. Observe that the claimed moments of A -1 are bounded: for all p < d, 

oo oo 

E [\- p ] < 1 + Ca Xj^k' 1 ^ = 1 + Ca k' 1 -^ log p (2 + k) < oo. 

k=l k=l 

For each fixed t > 1 + a\V\, the probability of encountering at least one ball 
centered in tV in which A < A|- t ] is 

(5.3) P [P k n tV ^ for some k > t] 

OO 

>1- 11 i 1 - atd \y\ k ~ 1 ~ d ) > 1 -exp(-a|V|). 

We now choose the constant a > to be large enough that the last quantity on 
the right, with \V\ = \Q\, is larger than the critical percolation probability p c (d) in 
dimension d. It suffices to take a := — 21og(l — p c (d)). 

Having constructed the random field A, we argue that homogenization fails for 
the equation (|5.2p by demonstrating that, almost surely, the solution u e blows up 
at the origin as s — > 0. 

Suppose that e < 1 and x £ and x/e 6 ^k>\i/e\Pki so that x/e is the 

center of a ball in which A < A[i/ e ] . By a simple rescaling argument we see that 

(5.4) inf u s {;u>) - mf u £ (-,w) > cXT^e 2 = ce|loge|. 

B(x,s/2) dB(x,e) /E 

In light of (I5.3[) , the probability of finding such a point x in each cube y + eQ is 
larger than p c {d). Since u £ is superharmonic, —Au e > in B±, we can apply the 
weak Harnack inequality to deduce that, in such a cube we have 

inf u £ > inf u £ + ce|loge| 

y+sQ/2 d(y+eQ) 

Observe that the collection {y + eQ : y £ (e/2)Z} of cubes tiles K d . The shape 
theorem in first passage percolation (c.f. Kesten (22j) implies that there exists 
a constant c\ > such that, almost surely and for sufficiently small e > 0, we 
can surround the origin by c|l/ej many disjoint walls of these good cubes of size e 
centered in (e/2)ZnQ. By a "wall" we mean a union of cubes of the form y+eQ with 
y G (e/2)Z such that the origin lies in a connected component of the complement. 
For each wall W, let dW = diW U d W where diW is the inner boundary (the 
boundary of the connected component of the compliment of W containing the 
origin) and d D W is the outer boundary (the boundary of the unbounded connected 
component of W). 

The construction above and another application of the weak Harnack inequality 
yield that, for each wall W, 

inf u e > inf u e + ce llogel . 

diW d a w 

Since there are c[l/ej such walls between the origin and dQ, we deduce that 

(5.5) u £ (0) > inf u e + (c|l/e|)e |loge| = c\ loge| -> +oo as e -> 0. 

dQ 

This verifies that both homogenization and regularity break down for (|5.2p . 

One wonders whether the failure of homogenization is simply due to the 1 on the 
right side of (|5.2[) . Our second example, which is a variation of the first, considers 
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an equation with zero right side in dimension d > 2 and with cllipticity having 
hnite pth moment for all p < d — 1. 

Write x = (x', x d ) e R d , with x' G K^ 1 and Q := (-1, l) d . Let A : R^ 1 -> (0, 1] 
be the random field constructed above, in one less dimension, and consider the 
Dirichlet problem 

(5.6) j-A(^)AV-<_=0 inQ, 

[ u £ = g on dQ, 

where A' is the Laplacian with respect to the first d—1 variables and < g < 1 is 
any continuous function on dQ' which is equal to 1 on {xd+i = ±1} and vanishes 
in dQ' ^{\xd+i\ < 1/2}. It is not difficult to show that the solution of (|5.6I) satisfies 
u e — > 1 locally uniformly in Q as e —> 0, which rules out homogenization to any 
effective elliptic problem. 

This confirms that homogenization may break down in dimensions larger than 
one for a linear operator even if the trace of its diffusion matrix is uniformly bounded 
below. In constructing the example, however, we destroyed the finite range depen- 
dence of the medium. We expect that this was necessary, since in the discrete 
setting Guo and Zeitouni [l9[ proved a quenched invariance principle for random 
walks in i.i.d random environments on the lattice, under the assumption that the 
random walker must move at each step (see also the interesting recent preprint of 
Berger and Deuschel Q which extends this to the non-strictly elliptic setting). We 
conjecture that, in the continuum setting, homogenization holds if we impose both 
a finite range dependence and require the diffusion matrix to have trace uniformly 
bounded away from zero. 



Appendix A. Elementary properties of the obstacle problem 

For the convenience of the reader, we briefly review (and sketch the proofs of) 
some well-known properties of the obstacle problem 



(A.l) 



mm {F(D 2 w,y),w} = inV, 
w = on dV. 



We have dropped the dependence of F on u) since the random environment plays no 
role here, and we furthermore assume that F is uniformly elliptic by the remarks 
preceding Theorem [TJ 

First of all, problem (|A.1[) has a unique solution w G C(V) which may be ex- 
pressed as the least nonnegative supersolution of F = in V: 

(A.2) w(x ; V) = inf {u{x) : u > in V and F(D 2 u, y) > in V) . 

This from the Perron method and the fact that the obstacle problem has a com- 
parison principle. Immediate from this expression is that w is a global subsolution: 

(A.3) F{D 2 w, y)>0 in V, 

as well as the monotonicity property: 

(A.4) V C W implies that w(- ; V) < w(- ; W) on V. 
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In order to use some regularity theory, we also need the fact that w satisfies 
(A.5) F{D 2 w,y) <k X{w=0} inV where k := sup(F(0, y))+. 

This is typically handled by considering an approximate equation with a penalty 
term (the Levy-Stampacchia penalization method) whose solutions satisfy (|A.5|) 
and showing that their uniform limit is w. We instead opt for a more natural and 
simpler proof by showing that w is also given by the formula 

(A.6) w(x ; V) — sup {u(x) : u < on dV and F(D 2 u, y) < kx{ u <o} in V] . 

It is clear that (|A.6|) implies (|A.5|) . To check the former, we let w denote the 
expression on the right side and observe that, since the zero function belongs to the 
admissible class by the definition of k, we have w > 0. The Perron method implies 
that w satisfies F(D 2 w,y) — in {w > 0}. Therefore w is a solution of (|4.1|) . and 
by uniqueness we deduce w = w. 
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